Theoretical study of even denominator fractions in graphene: Fermi sea versus paired 
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The physics of the state at even denominator fractional filhngs of Landau levels depends on the 
Coulomb pseudopotentials, and produces, in different GaAs Landau levels, a composite fermion 
Fermi sea, a stripe phase, or, possibly, a paired composite fermion state. We consider here even 
denominator fractions in graphene, which has different pseudopotentials as well as a possible four 
fold degeneracy of each Landau level. We test various composite fermion Fermi sea wave functions 
(fully polarized, SU(2) singlet, SU(4) singlet) as well as the paired composite fermion states in the 
n = and n = 1 Landau levels and predict that (i) the paired states are not favorable, (ii) CF 
Fermi seas occur in both Landau levels, and (iii) an SU(4) singlet composite fermion Fermi sea is 
stabilized in the appropriate limit. The results from detailed microscopic calculations are generally 
consistent with the predictions of the mean field model of composite fermions. 
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I. INTRODUCTION 

Although the fractional quantum Hall effect^ (FQHE) 
has not yet been observed in graphene, several papers 
have studied this possibility theoretically both in the 
SU(2) limit^ (when the Zeeman splitting is not small but 
the valleys are degenerate) and in the SU(4) limii^- (when 
both spins and valleys are degenerate). These studies 
show that, as for GaAs, the FQHE states are well de- 
scribed by the composite fermion (CF) theory^ and occur 
at the filling factors 



2pm ±1 



(1) 



where the partial filling of electrons (correspond- 

ing to total filling of — in — 2 + i/(")) or holes^ (at 
v = 4n + 2 — i/(")) in the graphene Landau level with 
index n, 2p is the CF vorticity, and m is the number of 
filled A levels (also known as CF Landau levels). There 
are differences, however. While FQHE in GaAs is much 
stronger in the lowest Landau level (LL), the FQHE in 
graphene is expected to be as strong in the |n| = 1 Lan- 
dau level as in the n = Ll^. More interestingly, many 
new incompressible CF states become possible because 
of the SU(4) symmetr}*2,. 

This work addresses the nature of the state at i/*^"'' = 
. If the model of weakly interacting composite fermions 
remains valid in the limit of ?7i ^ oo, then we expect a 
Fermi sea of composite fermions. In GaAs, the fully spin 
polarized Fermi sea of composite fermions has been ex- 
tensively studied^ and confirmed^ at ly ^ 1/2, and good 
evidence exists for a spin-singlet CF Fermi sea (CFFS) 
in the limit of vanishing Zeeman energy^'^. At v = 5/2 
in the second (n — 1) Landau level, it is currently be- 
lieved, although not confirmed, that the residual inter- 
actions between composite fermions produce a p-wave 
paired state of composite fermions, described by a Pfaf- 
fian wave functionl*^. In still higher Landau levels an 



anisotropic stripe phase is believed to occur. 

CF Fermi sea is an obvious candidate at half fillings in 
graphene, although it will have a richer structure asso- 
ciated with it. In the SU(4) symmetric limit, the mean 
field model of composite fermions predicts an SU(4) sin- 
glet CF Fermi sea, which has no analog in GaAs. The p- 
wave paired state of composite fermions is also a promis- 
ing candidate, especially at i^*^^) = 1/2 in the n = 1 LL, 
and it is interesting to ask if the graphene Coulomb ma- 
trix elements can make it more stable than the standard 
GaAs Coulomb matrix elements. For completeness, we 
also consider a so-called hollow-core state-'^^ describing 
the spin-singlet pairing of composite fermions, and, as in 
GaAs^«, find it not to be relevant. We note that our n — 
Landau level results below, as well as in Ref. [1, also ap- 
ply to the CF physics in valley degenerate semiconductor 
systemai^. 



II. MODEL 

The low-energy states of graphene are described 
in the continuum approximation by a massless Dirac 
Hamiltonian^^ 







(a-n)T 



AP, + gfiBB ■ S, (2) 



that acts on a 4-spinor Hilbert space. Here S denotes 
the spin and P the pseudospin associated with the valley 
degree of freedom, vp ~ 10® m/s is the Fermi velocity, 
H = p + |j4, and A is the on-site energy difference be- 
tween the two sublattices. The single particle spectrum 
of ifsr. is 



sgn(n) 
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Ap + g^isBs, (3) 



where s,p £ {57^5} s^ie the eigenvalues of Sz and Pz, 
respectively, and n is the Landau level index. In the 
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limit f/ ^ 0, A ^ each Landau level is 4-fold degener- 
ate, giving rise to an SU(4) internal symmetry. We con- 
sider below only the SU(4) symmetric part of the Hamil- 
tonian explicitly; from these results, the energy of any 
given wave function in the presence of certain kinds of 
symmetry breaking terms (for example, the Zeeman cou- 
pling) can be obtained straightforwardly, and level cross- 
ing transitions as a function of g and A can be obtained. 
The conditions for SU(4) symmetry have been discussed 
in Refs. H and[23l. 

Because we are interested in bulk properties, we will 
use the spherical geometry, in which electrons move on 
the surface of a sphere and a radial magnetic field is 
produced by a magnetic monopole of strength Q at the 
center. ■'^^'^^ Here 2Q(f)Q is the magnetic flux through the 
surface of the sphere; 4>o = hc/e, and 2Q is an integer 
according to Dirac's quantization condition. 

The interelectron interaction is conveniently 
parametrized in terms of pseudopotentialsi^ Vm,, 
where Vm is the energy of two electrons in relative 
angular momentum state m. The problem of interacting 
electrons in the n-th LL of graphene can be mapped into 
a problem of electrons in the n = LL with an effective 
interaction that has pseudopotential o^'^^ 



y(n)gr. _ 
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where the form factor Fn is 



Fo(fc) = 1, Fn{k) 



For an evaluation of the energies of various variational 
wave functions by the Monte Carlo method, we need the 
real-space interaction. In the n = LL this interaction is 
simply V{r) = l/r, where r is taken as the chord distance 
in the spherical geometry. In other Landau levels we use 
an effective real-space interaction in the lowest Landau 
level that reproduces the higher Landau level pseudopo- 
tentials in Eq. (|4]). We determine such an effective real 
space interaction in the planar geometry, and use it on 
the sphere. This procedure is exact in the thermody- 
namic limit, and it is usually reasonable also for finite 
systems. Following Ref. i3|, in the |n| = 1 LL we use the 
form 
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The coefficients {ci} are given in Ref. S We will assume 
parameters such that the finite thickness of the 2DEG 
and Landau level mixing have negligible effect. 

To build composite fermion trial wave functions, 
we will use the following consequence of Fock's cyclic 
condition^. The orbital part of one member of the SU(n), 
namely the highest weight state, can be constructed as 



\2p 



(7) 



where $s's are Slater determinants such that any state 
(n, m) in $s is also filled in ^s-i (conversely, if (n, m) is 
empty in $s, then it is also empty in $s+i); "Plll is the 
projection into the lowest {n = 0) Landau level^^; and 
the last factor, the Jastrow factor, attaches 2p vortices 
to each fermion. Here Ui — cos (6*2/2) e^"^*/^, and Vi — 
sin {Oi/2) e*"^'/^. The complete wave function is 



$'({r;})=^ $({,-}) n n 



(8) 
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where {a*} is a basis of the (n-dimensional) fundamental 
representation of SU(n), Mt is the number of particles 
in the a* state, mini = l,maxi = Afi,min2 = Mi -t- 
1, max2 = Ml + M2, . . . , and A is the antisymmetrizer. 

We define the CF Fermi sea as the thermodynamic 
limit of an integral number of filled Landau levels at an 
effective monopole strength q = for composite fermions. 
Clearly, if $1, . . . , $„ are identical, then Eq. ([7]) yields a 
legitimate trial wave function. We will label this state 
"CFFS . . . , ^]." As the effective monopole strength 
of composite fermions q is related to the real monopole 
strength Q as 

Q = q+p{N-l), (9) 
the filling factor is, assuming q = 0(1), 
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N^oo2Q + l Af^oo 2p(iV - 1) + 1 2p' 

The Pfaffian wave function^'^ , which is one of the can- 
didates for the FQHE state at = | in GaAs samplesi^, 
has the form 



^l/2p - 



UtVj - ViUj 
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on the sphere. By assumption, the Pfaffian wave function 
uses one spin band only. We also consider the hollow-core 



stati 



^hollow-core ^ (^^^^. ) "Q^^^^^. _ ^^.^^)2p^ (^2) 



i<j 



where Mij = (uiUj+jv/2~"i+JV/2^'j)^^ is t ^ t nratrix. 
This state is a spin singlet in the system with SU(2) 
symmetry; its symmetry becomes SU(2)xSU(2) in the 
SU(4) symmetric limit. Because of the last factor in Eqs. 
PT|) and (fT^ . which converts electrons into composite 
fermions, these wave functions describe paired states of 
composite fermions. 



III. RESULTS AND CONCLUSIONS 

We have studied CF Fermi sea states containing as 
many as 256 composite fermions (64 particles per Landau 
band), and our principal results^^ are given in Fig. [1] and 
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Table n These pertain to fillings 1^'-°'^ = \/2{v = ±3/2); 
1/(0) ^ 1/4 ^ ±7/4); zy(i) = 1/2 (i/ = ±5/2, ±11/2); 
1/(1) = 1/4 (i/ = ±9/4, ±23/4). (In relating i/(") to v, we 
have included the possibility of forming the state from 
either electrons or holes in a given Landau level.) When 
the spin or valley degeneracy is broken, the above study 
applies to many other half integral states also. To ob- 
tain the energy of the CF Fermi sea, we consider finite 
systems at _B* = and extrapolate the energy to the 
thermodynamic limit. The energies at v^^^ = 1/2 have a 
complicated dependence on l/iV, which makes extrapola- 
tion to the thermodynamic limit difficult. The following 
conclusions can be drawn. 

(i) For all fractions shown in Fig. [1] the hollow-core 
state has a very high energy and is therefore not relevant. 

(ii) The Pfaffian wave function also has higher energy 
than all of the CF Fermi sea states for all filling factors 
studied. In particular, it has higher energy than the fully 
polarized CF Fermi sea {[N]) in the n = 1 LL, in con- 
trast to GaAs where the fully polarized CF Fermi sea has 
higher energy^. We therefore conclude that the Pfaffian 
state is not stabilized in either n = or \n\ = 1 Landau 
level in graphene. Interestingly, for the fully polarized 
state, the overlaps given in Table |IT] indicate the Pfaf- 
fian wave function is actually a better representation of 
the exact Coulomb ground state at v^^^ — 1/2 in the 
n = 1 LL of graphene than it is of the 5/2 state in GaAs 
(for the latter, the overlaps are 0.87 and 0.84 for 8 and 
10 particles, respectively^^); nonetheless, energetic con- 
siderations rule out the Pfafhan state at i/^^^ = 1/2 in 
graphene. 

(iii) The overlaps given in Table HIl show that the Pfaf- 
fian is significantly worse at v^^^ = 1/2, indicating that 
it is not stabilized in the |n| = 2 LL of graphene either. 

(iv) We have considered CF Fermi sea wave functions 
of four distinct symmetries, ranging from SU(4) singlet 
to fully polarized. All of these have lower energies than 
either the Pfaffian or the hollow-core state. Without any 
symmetry breaking term, the SU(4) singlet CF Fermi sea 
has the lowest energy at v^^'' = 1/2, as expected from 
the model of non- interacting composite fermions. When 
the Zeeman and the pseudo-Zeeman energies are turned 
on, we expect a "partially-polarized" CF Fermi sea, and 
eventually a fully polarized CF Fermi sea. 

(v) The CF Fermi sea is also favored for = 1/2 
and i/(") = 1/4, but the energy differences between the 
various CF Fermi sea states are very small, less than the 
statistical error in our Monte Carlo evaluations. 

Other authors^i have considered a CF Fermi sea state 
at 1/ = 0, where the fourfold degenerate ti = LL is 
half full. Here, the electron (or hole) density in the 
ji — Landau level is p = 2\B\/(j)Q, which, upon com- 
posite fermionization of all electrons, gives an effective 



field oi B* = \B\ — 2(j)Qp = —3\B\ for composite fermions, 
which should be contrasted with B* — at 1/^"^ = 1/2. 
Khveschenko^i considers a state in which each of the 
four degenerate Landau bands is half filled forming a CF 
Fermi sea; the flux attachment docs not introduce corre- 



State 


£(.(") = i) 


£(.(«) = i) 


EiuW = i) 


CFFS [iV] 


-0.4651(1) 


-0.36014(4) 


n.a. 


CFFS [f,f] 


-0.46924(7) 


-0.35955(3) 


-0.3714(3) 


CFFS [f,f,f] 


-0.4732(1) 


-0.36019(6) 


-0.3720(2) 


CFFS [f,f,f,f] 


-0.47541(8) 


-0.36046(6) 


-0.3719(3) 


Pfaffian 


-0.45708(6) 


-0.35614(2) 


-0.3667(2) 


hollow-core 


-0.3141(3) 


-0.34932(3) 


-0.3564(2) 



TABLE I: The thermodynamic limit of the energy per parti- 
cle, in units of e^/el, for various CF Fermi sea (CFFS) states 
as well as the Pfaffian and the hollow-core wave functions at 
= i and i for |n| < 1. The notation for the CFFS state 
is explained in the text. 



N 


\n\ = 1 


H =2 


8 


0.902 


0.718 


10 


0.894 


0.486 



TABLE II: Overlap between the Pfaffian wave function *i/2 
and the exact ground state (the latter obtained assuming 
full spin and pseudospin polarization) at i/("' = 1/2 in the 
1 71 1 — 1,2 Landau levels of graphene in the spherical geome- 
try. (Note that the Pfaffian at A'^ = 6 and 12 occurs at the 
same flux values as v — 2/5 and 3/7, while A = 14 is beyond 
our computational ability.) 



lations between different bands in this approach. Finally, 
we comment on some of the approximations made in the 
model considered above. We have neglected LL mixing in 
our calculation; given that the energy difference between 
the CFFS and the Pfaffian states is fairly large ( 3-5%), 
we believe that LL mixing will not cause a phase tran- 
sition into a Pfaffian ground state, which is known to 
become worse with LL mixing^^. 
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FIG. 1: (Color online) Energy per particle, in units of e^/el, for several wave functions (four CF- Fermi-sea states, the PfafBan 
wave function, and the hollow core wave functions) at i/'"' = i and i in the n = Landau level (top) and in the |n| = 1 Landau 
level (bottom). Statistical error from Monte Carlo sampling is also shown. Here e is the background dielectric constant and I 
is the magnetic length. Extrapolation to the thermodynamic limit is given, wherever possible; the thermodynamic energies are 
quoted in Table |I1 
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try, the total energy is given by 



for all trial wave functions <l/t, where R is the radius of 
the sphere. The energy differences between various states 
do not depend on the details of the background subtrac- 
tion. Note that the graphene sheet in itself is not neutral. 



and part of the neutralizing charge resides on the back- 
gate in the experimental setup; the capacitive energy of 
the graphene sheet-backgate system is not taken into ac- 
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change their energy ordering. 
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